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ABSTRACT 

We report radio interferometric observations of the 12 C 16 1.3 mm J = 2 — 1 emission line in the 
circumstellar envelope of the M supergiant a Ori and have detected and separated both the SI and S2 
flow components for the first time. Observations were made with the Combined Array for Research 
in Millimeter-wave Astronomy (CARMA) interferometer in the C, D, and E antenna configurations. 
We obtain good u — v coverage (5-280 kA) by combining data from all three configurations allowing 
us to trace spatial scales as small as 0".9 over a 32" field of view. The high spectral and spatial 
resolution C configuration line profile shows that the inner SI flow has slightly asymmetric outflow 
velocities ranging from — 9.0 km s" 1 to +10.6 km s -1 with respect to the stellar rest frame. We find 
little evidence for the outer S2 flow in this configuration because the majority of this emission has 
been spatially-filtered (resolved out) by the array. We also report a SOFIA-GREAT CO( J = 12 - 11) 
emission line profile which we associate with this inner higher excitation SI flow. The outer S2 flow 
appears in the D and E configuration maps and its outflow velocity is found to be in good agreement 
with high resolution optical spectroscopy of K I obtained at the McDonald Observatory. We image 
both SI and S2 in the multi-configuration maps and see a gradual change in the angular size of the 
emission in the high absolute velocity maps. We assign an outer radius of 4" to SI and propose that 
S2 extends beyond CARMA's field of view (32" at 1.3 mm) out to a radius of 17" which is larger than 
recent single-dish observations have indicated. When azimuthally averaged, the intensity fall-off for 
both flows is found to be proportional to R , where R is the projected radius, indicating optically 
thin winds with p cx R~ 2 . 

Subject headings: circumstellar matter — radio lines: stars — stars: individual: (a Ori) — stars: 
late-type — stars: massive — supergiants 



1. INTRODUCTION 

The circumstellar envelope (CSE) of the M2 lab su- 
pergiant Betclgeuse (a Orionis) is a proving ground for 
ideas and theories of mass-loss from oxygen-rich M su- 
pergiants. Currently it is l osing mass at a respectable 
rate ~3x lCT 6 vr' 1 (iGlasseold fc Hugginsl [19861: 
iHuggins et al.1 119941 : lHarper et al.ll200lh . as it has been 
over the past ~ 1000 yr. Most of the opti cally thin sil- 
icate dust lies beyond ~ 46 stellar radii ([Danchi et al.l 
119941) and dust is, therefore, unlikely to be responsible 
for the bulk mass-loss. This raises the important point 
that if the mass-loss from Betclgeuse is not a result of 
dust then perhaps the same mechanisms that are respon- 
sible might also be active in the more dusty later M-type 
s upergiants. 

lHartmann fc Avrettl (|1984f ) constructed a WKB 
Alfven wave-driven model for Betelgeuse's chromosphere 
and wind that reproduced reasonably well the observed 
chromospheric emission fluxes and mass-loss rate. How- 
ever, multi-wavelength centimeter-continuum radio ob- 
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servations made with the Very Large Array (VLA) by 
iLim et aLI (|1998f ) revealed that the atmosphere was much 
cooler than the A lfven wave-driven model predicted 
(|Harper et al.l 1200 ID . Although new nonlinear Alfven 
wave models have been computed by lAirapetian et al.l 
(|2000[) . there are currently no theoretical models that 
make specific predictions for both the dynamic and ther- 
modynamic state of the mass outflow. Radiation pres- 
sure on atoms and molecules is another potential con- 
tributing candidate as a mass-loss mechanism and so 
spatial and dynamical studies of molecules are a fruit- 
ful line of investigation, especially in relation to eventual 
formation of dust. Such studies also allow us to calculate 
the time scales on which certain mass-loss episodes have 
occurred, and these can then be compared to the time 
scales of potential mass-loss initiators such as convection 
or magnetic dynamo cycles. 

The study of CO molecules in the CSE of Betelgeuse 
began with the detection of 4.6 /im r o- vibrational absorp - 
tion lines of 12 C 16 and 13 C 16 bv lBernat et ill 0fa 
who identified two absorption features, implying two dis- 
tinct structures within the overall outflow. One compo- 
nent, known as SI, has a Doppler shift of 9km s _1 toward 
us with Texc — 200K, Vturb — 4kms~ 1 and 7Vi2 C i6 = 
4.7 x 10 17 cm~ 2 . The second faster component, known 
as S2, has a Doppler shift of 16 km s _1 towards us with 
T cxc ~ 70 K, Wturb — lkms -1 and Ni^qwq = 1.2 x 
10 16 cm -2 . The SI feature with its higher column density 
was well known from atomic absorption line studies (e.g., 
lWevmannlll962D and both features had been detected in 
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Table 1 



CARMA Observations 


Observation 


Configuration 


Time on Source 


Flux 


Phase 


Image Cube a 


Date 




(hr) 


Calibrator 


Calibrators 


Dynamic Range b 


2007 Jun 18 


D 


0.9 


0530+135 


0530+135, 0532+075 


22.8 


2007 Jun 21 


D 


3.0 


0530+135 


0530+135, 0532+075 


22.7 


2007 Jun 24 


D 


2.1 


0530+135 


0530+135, 0532+075 


26.1 


2007 Jun 25 


D 


2.4 


0530+135 


0530+135, 0532+075 


30.2 


2009 Jul 7 


E 


3.2 


3C120 


3C120, 0532+075 


30.1 


2009 Nov 5 


C 


1.2 


3C120 


3C120, 0532+075 


17.3 


2009 Nov 09 


C 


3.0 


3C120 


3C120, 0532+075 


27.2 


2009 Nov 15 


c 


1.0 


3C120 


3C120, 0532+075 


17.8 


2009 Nov 16 


c 


3.2 


3C120 


3C120, 0532+075 


32.0 


All 


c 


8.4 






43.8 


All 


D 


8.4 






31.9 


All 


Multi-configuration 


20.0 






52.3 



Low spectral resolution (i.e., channel width of 1.3 km s ). 

The peak emission of the image cube divided by the root mean square of the residual image. 



high spectral resolution ato mic Na and K absorption pro- 
files (|Goldberg et al.lll975h . 12 C 16 was subsequently 
detected at 230 GHz in the J = 2 — 1 rotational emis- 
sion line by IKnapp et al.1 (I1980D. although a searc h for 
SiO(J = 2 - 1) bv lLambert fc Vanden Boutl ([1971 had 
been unsuccessful. The we aker 12 C 16 Q(J = 1 — ) line 
was tentatively detected bv IKnapp fc Morrisl (|1985| ) with 
a 7m dish which had a half-power beamwidth (HPBW) 
o f 100". 

iHugginsI (|1987|) presented a higher signal to noise 
12 C 16 0( J =2-1) observation of Betelgeuse's CSE with 
an HPBW of 32" and found some evidence for an S2 
radius of about 16" by comparing the (2 — 1)/(1 — 0) in- 
tensities. However, a 30 m Institut de Radioastronomie 
Millimetrique (IRAM) J = 2 — 1 line profile was later 
presented by iHuggins et al.1 (|1994[ ) and looked remark- 
ably similar, even though it was observed with a smaller 
12" HPBW. The profile did not show the horned wing 
signature expected if it had been resolved apparently in 
conflict with the previous S2 radius estimate. 

Here we present the results of an interferometric study 
of the rotational 12 C 16 0(J = 2 — 1) emission line made 
using three Combined Array for Research in Millimeter- 
wave Astronomy (CARMA) configurations with HPBWs 
of 0".9 (C), 2".l (D), and 4".4 (E) designed to explore 
the SI and S2 flows at these spatial scales. Prelim- 
inary results of the D configuration observations have 
been presented in lHarper et al.l (|2009[ ) as part of a multi- 
wavelength study of CO surrounding a Ori. We also 
present a supporting SOFIA CO(J = 12 — 11) line pro- 
file in addition to high spectral resolution observations 
of the K 1 7699 A line. In Section 2 the observations and 
data reduction techniques are discussed and in Section 3 
the results of the spectra and image maps arc presented. 
Discussions and conclusions are presented in Section 4 
and Section 5, respectively. 

2. OBSERVATIONS AND DATA REDUCTION 

The millimeter observations were made with th e 15 el- 
ement CARMA interferometer (jScott et al.ll2004| ) which 
is located at Cedar Flat in eastern California at an eleva- 
tion of 2200 m. The array consists of nine 6.1m antennas 
and six 10.4m antennas formerly from the Berkeley Illi- 
nois Maryland Association (BIMA) and the Owens Val- 
ley Radio Observatory (OVRO) arrays respectively Ta- 
ble [1] summarizes the various observations which span 



the period 2007 May— 2009 November. The observa- 
tions were carried out in the C, D, and E configura- 
tions and consist of on-source profiles of the 12 C 16 0( J = 
2 — 1) line which has a rest frequency of 230.538 GHz 
(1.3mm). The baseline length spans over 26— 370m (C 
array), 11 — 148m (D array), and 8.5— 66m (E array) pro- 
viding HPBWs of 0".9, 2".l, and 4" .4 respectively at 
1.3mm. The HPBW of the individual 10.4m antennas 
is ~ 32" at the observed frequency. 

The CARMA correlator takes measurements in three 
separate bands, each having an upper and lower side- 
band. One band was set to the low-resolution 468 MHz 
bandwidth mode (15 channels of 31.25 MHz each) to ob- 
serve continuum emission and was centered on the line. 
The other two bands were configured with 62 MHz and 
31 MHz bandwidth across 63 channels (with a resolution 
of 1.3km s _1 and 0.65km s~ , respectively) and were also 
centered on the line. The line was measured in the upper 
sideband in the C and E array and in the lower sideband 
in the D array. 

Bandpass and phase calibration were performed using 
3C120 and 0530+135. 0532+075 was used as a secondary 
phase calibrator to determine the quality of the phase 
transfer from the primary phase calibrator. The observ- 
ing sequence was to integrate on the primary phase cal- 
ibrator for ~ 2.5 minutes, the target for ~ 18 minutes, 
and the secondary phase calibrator for ~ 2.5 minutes. 
The cycle was repeated for each track which lasted be- 
tween 1.5 hr and 5 hr. Absolute flux calibration was 
carried out with 0530+135 and 3C120 using the continu- 
ously updated CARMA flux catalog to obtain their flux 
values at each observation. 

The raw data were smoothed by a Hanning filter within 
MIRIAEH and then exported into FITS format so that it 
could be analyzed with the CASA0 data reduction pack- 
age. All calibration and imaging was carried out within 
CASA. The image cubes were multi-scale CLEANed 
down to the 3er threshold using natural weighting and 
were corrected for primary beam at tenuation. The mul- 
tiscale algorithm (jRich et all 120081 ) within CASA was 
set to four unique scales; the largest corresponding to 

7 Multichannel Imag e Reconstruction, Im age Analy sis and Di s- 
play, http : //www. atnf . csiro . au/computing/ sof tware/miriad/ 

8 Common Astronomy Software Applications, 
|http://c 



/casa.nrao 



edu/| 
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the largest structures visible in individual channel maps. 
Each scale was approximately set to three times smaller 
than the preceding scale. 

Each of the three CARMA configurations sample a dif- 
ferent range of spatial frequencies; the range of which is 
dependent upon the maximum and minimum baselines 
(bmax and 6 m i n ) of each configuration. The sources we 
are observing are extended and therefore it is necessary 
to consider the response of each CARMA configuration 
to this extended emission. For any array configuration, 
emission with angular scale s of ~ A /b m in or g reater is 
not reproduced in the maps (jTavlor et al .11 19991 ) and this 
scale is often used as a guide for the resolving out scale 
or maximum scale of an array configuration. To obtain 
a more robust estimate of the largest angular scale that 
can be accurately imaged in the high spatial resolution C 
configuration maps we computed the visibilities of an ex- 
tended emission feature (whose spatial extent was set to 
that of the primary beam) using CASAs simulation tool, 
simdata. This tool then produced a CLEANed image of 
these visibilities from which we calculated the resolving 
out scale to be ~ 6" (i.e., 0.6A/& m ; n ). Ultimately, com- 
bining the data from these three configurations allows 
the missing short spacings from the extended C configu- 
ration to be recovered while maintaining its high spatial 
resolution. 

3. RESULTS 

Bctclgcusc is a semi-regular variable and its radial ve- 
locity exhibits variability on time s cales ranging from 
short 1.5 year periods as suggested bv lStebbins fc Huffed 
(fl93l to Ion ger 5.8 year periods (Jones, 1928). Its ra- 
dial velocity amplitudes are also known to vary by at 
least ±3 km s _1 (jSmith et al.lll989l ) making it difficult to 
determine a precise value for the stellar center-of-mass 
radial velocity. In this study we adopt a heliocentric ra- 
dial velocit y of +20.7 km s -1 fa sr = 4.8 km s _1 ); a value 
adopted bv lHarper et al.l (120081) and i s base d on the mean 
values of Jones (1928) and iSanfordl (|1933[ ). All spectra 
are plotted with respect to the stellar center-of-mass rest 
frame. 

3.1. CO(J =2-1) Spectra 

The spectrum for each individual configuration image 
cube (which are composed of all the appropriate config- 
uration tracks listed in Table [1} along with the multi- 
configuration image cube can be used to obtain infor- 
mation on the kinematics of the SI and S2 flows. The 
spectra corresponding to the C, D, and E configuration 
image cubes are plotted in Figure Q] for both the high 
(0.65 km s" 1 bin" 1 ) and low (1.3 km s" 1 bin -1 ) spectral 
resolution data and were obtained by integrating all emis- 
sion within a circular area of radius 5" centered on the 
source. The high and low spectral resolution modes al- 
low two independent sets of spectra to be measured for 
each observation and thus provide a good check on the 
data quality. The high resolution spectra (channel width 
= 0.65 km s~ ) give the best measure of SI and S2 kine- 
matics and therefore all outflow velocities are derived 
from these spectra. 

The E configuration image cube spectrum has a total 
line width of 29.2 kms -1 and the low spectral resolu- 
tion profile contains a steep blue wing emission feature 



1 1 1 1 1 1 1 

1.3 km s" 1 
0.65 km s" 1 
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Figure 1. Spectra integrated over a radius of 5" for each array 
configuration image cube. The blueshifted emission component be- 
tween — 16.0km s — 1 and —10.0km s _1 is almost resolved out in the 
C configuration image cube spectrum. The red and blue lines cor- 
respond to the high and low spectral resolution data respectively. 

between —16.0 km s -1 and —11 km s _1 and a more flat- 
topped feature between — 10.3kms _1 and — 13.2kms . 
This steep emission wing shows that the turbulence in 
the flow is less than or equal to the velocity bin size. 
The blue wing in the high resolution profile matches the 
lower resolution profile well but the remainder of the pro- 
file looks more complex than the fiat-topped feature seen 
in the lower resolution profile. The profile shape of the 
CO (J = 2 — 1) line has been well documented by pre- 
vious single-d ish observations (e.g., IKnapp et al.l 119801 
IHugginsj|1987t) and, out of our three individual configu- 
ration spectra, we expect the most compact E configura- 
tion spectra to resemble these single-dish measurements 
the closest due to its better sampling of the inner u — v 
plane and consequent sensitivity to extended structures. 
This indeed turns out to be the case when we compare 
our three individual configuration spectra to those pre- 
vious single-dish profiles. The blue wing emission fea- 
ture appears again in the D configuration spectrum at 
the same velocities as those in the E configuration spec- 
trum but the remainder of the profile appears quite dif- 
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ferent. Between — 10.3 km s and +13.2 km s the D 
configuration spectrum is dominated by a blue wing at 
~ — 10.0 kms" 1 , a red wing at ~ +13.0 kms -1 and an 
emission feature at ~ km s" 1 . 

The line profile has a much lower flux in the high spa- 
tial resolution C configuration spectrum due to its lack 
of sensitivity to extended structure. The blueshifted 
emission feature located between —16.0 kms" 1 and 
— ll.Okms -1 in the E and D configuration spectra is 
almost completely spatially filtered by the extended C 
configuration. This component of the line has previously 
been associated with the outer S2 flow dHugginsI 119871) 
and as the majority of it has been spatially filtered by 
our C configuration we expect even less contribution from 
the S2 flow at lower absolute velocities still. For the red- 
shifted line emission we again expect the majority of the 
S2 contribution to be spatially filtered, so we conclude 
that the majority of the emission in the C configura- 
tion spectrum emanates from the inner SI flow. The 
spectrum is double peaked with the blue and redshiftcd 
wings extending to — 9.0 kms -1 and 4-10.6 km s" 1 re- 
spectively, and we define these as the outflow velocities 
of SI. As discussed in Section 2, the C configuration has 
a resolving out scale of ~ 6" at 1.3 mm and so is not sen- 
sitive to angular scales larger than this. If the emission 
between — 9.0 kms" 1 and +10.6 kms" 1 in the C config- 
uration spectrum appeared as a flat-topped profile then 
we could conclude that the SI flow lies within a radius 
of 3" from the star. Clearly however, the lower absolute 
velocity components of this profile have been spatially 
filtered so we conclude that the radial extent of the SI 
from the star is greater than 3". If we assume that the 

51 flow would produce a top-hat line profile were it not 
for the resolving out effects of the interferometer, then 
its integrated line flux is 3.1 x 10 -19 Wm" 2 . 

To obtain the most robust value for the S2 outflow 
velocities we examine the high spectral resolution multi- 
configuration image cube spectrum which is composed of 
all tracks from all three configurations. It is worth stress- 
ing that by analyzing the multi-configuration image cube 
we make the assumption that the physical properties of 
all three components (i.e., a Ori, SI, and S2) have not 
changed over the total observation period (i.e., ~ 2.5 yr). 
The profile is found to have a total line width of 28.6 
± 0.7 kms" 1 , which is in close agreement with previous 
single-dish observations of the line where values of 30.6 ± 
2.5km s _1 a nd 28.6km s -1 w ere reported bv lKnapp et al.l 
(|1980l) and IHugginsI (|1987l ). respectively. The centroid 
velocity of the spectrum is -1.1 ± 0.7 kms" 1 (uj sr = 3.7 
± 0.7 km s~ ) whi ch is again in close agreement with 
iKnapp et all (|1980| ) and IHugginsI (|1987f ) values of u lsr = 
3.0 ± 2.5 kms" 1 and vi sr = 3.7 ± 0.4 kms -1 respec- 
tively. The integrated line flux is 1.5 xlO" 18 Wm" 2 of 
which approximately 20% emanates from the SI flow. 

The outflow velocities of S2 are —15.4 kms -1 and 
4-13.2 kms -1 which, like the SI flow, are slightly asym- 
metric but in the opposite sense. Note that the SI and 

52 outflow velocities are dependent on the adopted ra- 
dial velocity of Betelgeuse. If, for instan ce, we instead 
adop t a radial velocity of 21.9 km s" 1 dFamaev et al.l 
120051 ) then the S2 outflow velocities become even more 
asymmetric (—16.6 and +12.0 kms" ) while the SI out- 
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Figure 2. Spectral profiles of the low spectral resolution multi- 
configuration image cube for circular extraction areas of radius 1", 
2", 4", 6", 8", and 10". The signal to noise of the line profile 
reduces at larger extraction areas as more noise is included from 
the outer regions of the channel maps. 

flow becomes less so (—10.2 and +9.4 kms" 1 ). Both SI 
and S2 therefore cannot have spherically symmetric out- 
flow velocities regardless of the adopted stellar radial ve- 
locity. Adopting a m ass of 18 M Q and a radius of 950 Rq 
(|Harper et al.l 120081 ) then the escape velocity for Betel- 
geuse is 85 kms" 1 which is much greater than the SI 
and S2 outflow velocities. This indicates that the major- 
ity of the stellar mass-loss mechanism's energy goes into 
lifting the CO molecules out of the gravitational poten- 
tial and not into their outflow velocities. These outflow 
velocities are greater than the adiabatic hydrogen sound 
speed, which, if we assume that the gas temperature is 
the same as the excitation temperature, are 1.7 kms" 1 
and 1 kms" 1 for SI and S2, respectively. 

The spectra in Figure [2] are taken from the low- 
resolution multi-configuration image cube using circular 
extraction areas ranging in radius from 1" to 10" and 
demonstrates how the line profile changes over these dif- 
ferent extraction areas. The most striking change in the 
line profile is the change in appearance of the extreme 
blue wing. At small extraction radii where we sample 
the most compact emission, the feature is weak in com- 
parison to the rest of the line but becomes more domi- 
nant as we begin to sample more of the extended emis- 
sion. This indicates that even the high absolute velocity 
components of the S2 flow have extended emission and 
this is why they are almost completely spatially filtered 
by CARMA's C configuration. The large reduction of 
flux at —11 kms -1 suggests that there is more material 
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moving toward the observer than at other lower absolute 
velocities indicating a non-isotropic (or non-spherical) S2 
flow. This suggests a more sheet-like (flatter) structure 
rather than a spherical cap. 

3.2. Multi- configuration Image Cube 

A subset of the blueshifted velocity channel maps of the 
low spectral resolution multi-configuration image cube 
is presented in Figure El The first channel map at 

— 17.9kms _1 shows just the compact unresolved con- 
tinuum emission with no extended emission present. Be- 
tween — 16.7km s _1 and —9.0km s _1 , we see evidence for 
the development of a classical shell signature for the S2 
flow. We first sample the highest velocity components 
where the emission is relatively compact (i.e., between 

— 16.7kms _ and — 12.9 km s~ ) and then sample lower 
radial velocity components where S2 becomes a faint ring 
(i.e., between — 11.6 km s -1 and — 9.0kms _1 ). At lower 
velocities again, these rings disappear into the noise of 
the maps and possibly extend out beyond the primary 
beam at zero velocity when the rings should have maxi- 
mum spatial extent. The emission from the channel maps 
between — 15.3 km s -1 and — 11.6 km s" 1 corresponds to 
all the emission in the extreme blue wing component of 
the multi-configuration image cube line profile discussed 
in Section 3.1. We can see in Figure [3] that all of this 
emission is greater than the C configuration resolving 
out scale therefore confirming that our C configuration 
line profile is mainly composed of SI emission. The shell 
formation signature of S2 is also apparent in the red- 
shifted velocity channel maps between +7.5 kms -1 and 
+13.8 km s" 1 but the emission appears weaker and the 
rings fainter therefore indicating that S2 is somewhat 
fragmented. 

The multi-configuration maps also show the central 
compact emission from the SI flow at velocities between 
— 10.3 km s -1 and +11.3 km s~ . This SI emission can 
be seen in the final two maps of Figure Has a central 
slightly elongated emission feature surrounded by the 
fainter rings of the S2 flow. In the maps where both SI 
and S2 are present the emission from SI appears brighter 
than the emission from S2. The spatial extent of the SI 
flow varies from channel map to channel map but ap- 
pears to be larger than the 2" value given bv lSmith et al.1 
()2009[) , who observed off-star wind scattered ro- vibration 
CO lines. 

An additional spatially unresolved source is detected in 
a number of the D configuration image cube maps (both 
high and low spectral resolutions) and ha s been previ- 
ously documented by lHarper et al.l (j2009f ). The com- 
ponent is present in only five continuous low-resolution 
channels between ~ — 4.0kms~ 1 and +2.4 km s" 1 and 
is located ~ 5" S — W of a Ori as shown in Figure HI 
Its peak emission lies at ~ kms -1 and here it approx- 
imately equals 60% of the source peak emission. The 
corresponding channel maps in the E configuration im- 
age cube show extended emission out to 8" in the same 
S — W direction. This second source does not appear in 
any of the C configuration channel maps probably due to 
the lower sensitivity resulting from the smaller HPBW 
(i.e., the flux is diluted). This discrete second source 
thus has the effect of adding extra emission to the cor- 
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-0.1 +0.1 +0.2 

Figure 3. Eight channel maps from the multi-configuration image 
cube (Av = 1.3km s ). The peak emission has been cut at 0.2 Jy 
beam - to emphasize the fainter emission. The color scale is linear 
and has been normalized to this maximum cutoff and minimum 
value of each channel. The contour levels are at —2a, 2a, Aa, and 
6a (la ~ 0.03 Jy beam -1 but varies per channel). Dashed green 
lines represent negative contours. 

responding multi-configuration image cube maps at the 
low velocities where it is present. 

3.3. Determination of the SI and S2 Radii 

The spatial extent of the SI and S2 flows around Bctel- 
geuse was not directly determi ned from eith e r the CO 
infrared absorption spectra of IBernat et al.l (I1979T) or 
previous CO single-dish radio observation s ()Knapp et al.l 
ll980tlHuggrnill987l ; lHuggins et alj|1994h . Our low spec- 
tral resolution multi-configuration image cube has suffi- 
cient spatial resolution and signal to noise to make direct 
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Figure 4. Integrated intensity image of the D configuration chan- 
nel maps that contain the discrete second source approximately 5" 
S — W of a Ori. Contours for the integrated intensity are la, 1.5a, 
la, and 3a (la = 1.3 Jy beam -1 km s -1 ). The size of the restoring 
beam is shown in white in the bottom left corner. 

estimates of the maximum radius of both flows. The 
outer S2 flow is not seen in the low absolute velocity 
channel maps where its spatial extent is maximum and 
either lies outside of the primary beam or is lost into the 
noise near the edge of the maps. We derive the maxi- 
mum outer scale of the S2 flow by looking at the spatial 
scales of the S2 flow in the higher absolute velocity maps 
where S2 is present. If we assume that S2 is spherically 
symmetric with an outer radius Rs2, and is undergoing 
steady expansion with velocity V<32, then we can estimate 
its radius in each velocity channel using the following re- 
lation: 



^*cha 



i? S 2sin 



^chan 



(1) 



where r c h an is the S2 radius in a channel at velocity u c han- 
We use Equation (1) to estimate the maximum pro- 
jected spatial extent of S2 which occurs at zero velocity 
An estimate of the S2 radius per channel (r c han) was 
found by creating annuli of increasing radius around the 
central emission in each relevant line channel map of the 
multi-configuration image cube, extracting all flux within 
each annulus and then plotting these fluxes against dis- 
tance from the star for each channel. The maximum of 
these resultant curves was then deemed to be the max- 
imum radius of S2 per channel. Figure [5] shows these 
data over-plotted with two model outflows which were 
created using Equation (1). The blueshifted data points 
were best fitted by a model outflow of maximum radius 
17" and outflow velocity 17 kms -1 , while the redshiftcd 
data points were best fitted by a model outflow of max- 
imum radius 16" and outflow velocity 14kms~ . It is 
worth mentioning that this estimate for the spatial ex- 
tent of S2 is only weakly dependent on our adopted radial 
velocity value for Betelgeuse and adopting a slightly dif- 
ferent value would simply alter S2s outflow velocities. As 
S2 is not present in our lowest absolute velocity map we 
are not able to report an estimate of its width. 

In the left column of Figure [6] we investigate the in- 
tensity distribution of CO emission as a function of pro- 
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Figure 5. Derived S2 radius as a function of velocity (red points) 
overplotted with two model outflows. The blueshifted model (left) 
corresponds to an outflow with a maximum radius of 17" and a 
velocity of 16.7 km s" 1 while the redshifted model (right) corre- 
sponds to an outflow with a maximum radius of 16" and a velocity 
of 13.8 kms" 1 . Note: the line profile is 1.9 1.9 kms" 1 wider in 
the low resolution image cube (Av = 1.3 kms -1 ) than in the high 
resolution image cube. 

jected radius, R, for both the SI and S2 flows. From 
our discussions in Section 13.11 we can assume that all 
line emission between —15.4 — > — 10. 3 km s -1 and +12.4 
— > + 13.8 kms -1 emanates solely from the S2 flow. Us- 
ing the low spectral resolution multi-configuration image 
cube we integrate the surface brightness over these chan- 
nels and find that the intensity fall off is proportional to 
i? -1 (Figure El top). To investigate the SI flow intensity 
distribution around a Ori we integrate the surface bright- 
ness over the channels between —9 — > +10.6 km s -1 . Al- 
though these channels contain emission from both SI and 
S2, most of the S2 emission here will have larger pro- 
jected radii and thus the majority of the inner emission 
should emanate from the SI flow. Between 0".5 and 4" 
from the star the intensity is again found to be propor- 
tional to R^ 1 (Figure [U bottom). Such an intensity dis- 
tribution is expected for an optically thin homogeneous 
constant velocity outflow with p oc 1/R 2 . Beyond 4" the 
intensity fall off is more rapid and is close to an R~ 2 dis- 
tribution which may mark the initiation of the current 
epoch of mass-loss. 

Insight can also be gained into how the intensity varies 
on different size scales by conducting analysis in the 
u — v plane and plotting the visibility amplitude of a Ori 
against u — v distance. The result of this is shown in the 
right column of Figure [6] where the same channels cor- 
responding to the SI and S2 flows have been used. The 
data arc azimuthally averaged, and have been binned to 
produce one data point per kA. The result for both the SI 
and S2 data is a steep drop-off in visibility amplitude over 
a relatively short u~v distance, signaling that the sources 
are well resolved. Both sets of visibility data agree with 
an intensity proportional to (a 2 + i? 2 ) -1//2 , where a is an 
inner spatial limit. This is be cause the Hanke l transform 



of this function is q 



-2naq 



(|Bracewellll2000l ) , where q is 



the u — v distance, and a vertically scaled version of this 
function is shown to match the visibility data very well in 
Figurc[6l As analysis in both the sky and u— v plane indi- 
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Figure 6. Left Column: surface brightness as a function of projected radius on sky, R (red line). The emission has been extracted 
from the low spectral resolution multi-configuration image cube and is integrated over the channels where SI is present (bottom) and over 
the channels where only S2 is present (top). Intensity proportional to R~ 1 and R~ 2 is also shown for comparison. Right Column: the 
corresponding visibility amplitude as a function of u — v distance (q) of both outflows can be modeled well by an R^ 1 fall off in intensity. 
The error bars in all plots represent the standard error of the mean. 



cate the intensity of both flows is proportional to _R _1 we 
conclude that when azimuthally averaged, both outflows 
are consistent with an opticall y thin and quasi-s teady 
flow whi ch is in agreement with iSmith et al.l (120091 ) (i.e., 
SI) and lPlez fc Lambertl ((2002T) (i.e., S2). 

An exact determination of the maximum spatial ex- 
tent of the SI flow is more difficult as we do not see 
the classical shell formation signature for it as we sample 
across velocities, like we do for S2. Instead its spatial 
extent varies over the channel maps with evidence of dis- 
crete clumps being present in many of these maps. At 
20% of maximum emission in the integrated intensity 
SI map (i.e., composed of all channels between — 9 — > 
+10. 6km s _1 ) the SI flow extends out to a mean distance 
of ~ 4" and is even more extended in the S — W direc- 
tion due to the presence of the second emission feature 
in the compact configuration data sets. The HPBW of 
0".9 is not sufficient to determine whether the SI flow 
is discrete or an extension of t he current wind phase 
seen i n ultraviolet spectra, e.g., iCarpenter fc Robinson! 
( 1997T) . and centimeter-radio contin uum interferometry 
(jLim et al.lfl998t lHarper et al.ll2001[ ). 



3.4. Continuum Flux Densities 

In Table 2 we show the derived continuum flux densi- 
ties for each of the three configuration image cubes and 
also the multi-configuration image cube. The high spec- 
tral resolution (Av = 0.65 km s -1 ) image cubes were just 
wide enough to image the CO line but were too nar- 
row to make accurate estimates of the continuum flux 
density. Therefore, all continuum flux density estimates 
are derived from the lower spectral resolution (Aw = 
1.3kms _1 ) image cubes from which we were able to take 
accurate measurements at both sides of the line. We 
fitted elliptical Gaussians to ~ 20 continuum channels 
using CASA's imfit routine allowing the flux and corre- 
sponding uncertainties to be calculated. The source was 
unresolved in most of these continuum channels. 

Betclgeuse is kn own to show bri ghtness variations at 
many wavelengths. I Goldberg! (|1984l ) reports a decrease of 
half a magnitude in visua l bright ness over a period of six 
years. iBookbinder et ail (|1987t ) found stochastic 30%- 
40% variations in flux density at 6 cm over timescales as 
short as 10 days to as long as 8 months (i.e., the ob- 
servational per iod). A mo r e com prehensive study was 
carried out by iDrake et al.l (|1992f) who observed Betel- 
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Table 2 

CARMA Continuum Fluxes at 230 GHz 



Configuration 


Restoring Beam 


Flux 


Uncertainty 




(" X ") 


(mjy) 


(mjy) 


C 


0.96 x 0.76 


234 


18 


D 


2.33 x 1.87 


389 


72 


E 


4.93 x 3.84 


278 


10 


Multi-configuration 


1.05 x 0.84 


289 


21 



geuse with the VLA at centimeter-wavelengths from 
1986 to 1990 and found stochastic variability of 22%, 
15%, and 21% at 6 cm, 3.6 cm, and 2 cm, respectively 
at a variety of different timescales down to less than 1 
month. The millimeter-continuum emission that we mea- 
sure arises mainly from electron— ion and electron— atom 
brcmsstrahlung and possibly dust emission, so it is not 
unreasonable to also expect variability at millimeter- 
wavelengths too. The D configuration data were ac- 
quired under adverse weather conditions and these data 
have the highest noise levels out of the three configura- 
tions. Its continuum emission measurement is approxi- 
mately 50% greater than the C and E configuration con- 
tinuum measurements which were also acquired approx- 
imately two years after the D configuration data. We 
believe that the continuum emission derived from the 
multi-configuration image cube is a reasonable estima- 
tion of the mean millimeter-continuum flux density over 
the two year period and is in reas onably good agr e ement 
with the 250 GHz flux density of lAltenhoff et all ([l994h 
who report a mean value of 351±25 mjy for 1986 — > 1989. 

4. DISCUSSION 
4.1. Previous CO Observations 

iBernat et al.1 (|1979D were the first to detect circum- 
stcllar absorption lines in CO by looking at the fun- 
damental ro-vibration lines at 4.6 fira. These infrared 
observations revealed two separate outflows around a 
Ori; an excited (T oxc = 200 K) SI flow with an expan- 
sion velocity of 9 kms -1 and a less excited (T exc = 70 
K) S2 flow moving with a faster expansion velocity of 
16kms _1 . iKnapp et all (fl980h were the first to detect 
emission in the CO(J = 2 — 1) line at 1.3 mm using 
the 10 m millimeter- wave telescope at OVRO but only 
detected one component expanding at 15 km s _1 . By an- 
alyzing the shape of the line pr ofile, they con c luded that 
the S2 radius of 55" derived by IBernat et al.l (|1979f ) was 
too large and that it lies at a radius of R < 10". Since 
the detection by Knapp, a number of observations at 1.3 
mm have been carried out with various beam sizes and all 
spectra look remarkably similar; that is the profile has a 
steep extreme blue shifted emission component with the 
remainder of the profile looking more flat-top ped and 
contai ning a number of less dominant spikes. iHugginsI 
<|1987D used their single-dish observations (HPBW ~ 32") 
of the CO(J = 2 — 1) line along with excitation and 
self-shielding models of CO to conclude that the SI flow 
makes little contribution to the final emission line. They 
also identify the extreme blue wing of the line with the 
S2 flow and predict that it may extend out to a radius of 
~ 16". Later however, iHuggins et al.l (|1994j ) compared 
their detected 609 fim 3 Pi — > 3 Pq fine structure line of C 
I with CO data obtained with the IRAM 30 m telescope 



(HPBWs ~ 12") and find that the expansion velocities in 
both lines are essentially the same. They conclude that 
the radial extent of C I is < 7" and both the CO and 
C I are formed in the inner envelope and roughly extend 
over the same area. 

The shape of our multi-configuration line profile for ex- 
traction areas of radii 6" or greater are in good agreement 
with previous hi gh signal to noise sin gle-dish CO( J = 2 — 
1) spectra (e.g. IHuggins et al.|[l994l Figure 1) although 
the emission spikes in our line profiles are more dominant. 
Our total line width of 28 .6 km s" 1 is in g ood a greement 
with IHugginsI (|1987l) and IHuggins etall (|1994T > who re- 
port line widths of 28.6 kms -1 and 30 kms -1 respec- 
tively. The extreme blue wing in both of these spectra 
is the dominant emission feature of the line and this is 
also true in our multi-configuration spectra at extrac- 
tion areas > 6". Using the IRAM 30 m te lescope, which 
has an HPBW of only 12" at 230 GHz, IHuggins et all 
(|1994p pr oduce a s i milar line profile shape to that pre- 
sented in IHugginsI (|1987l ) who have a larger HPBW of 
~30". From this, one would expect that the major- 
ity of the blue wing emission is compact. Our multi- 
configuration line profiles suggest otherwise however, and 
show a continuous increase in the blue wing emission as 
we take larger extraction regions out to 10". The multi- 
configuration maps also show a faint ring structure form- 
ing at ~ —11.6km s _1 and expanding further out in lower 
absolute velocity channel maps. This ring emission is 
fainter than the higher velocity compact emission so we 
see a drop in flux density in our spectra at the point 
where these rings form. Therefore, the steepness of the 
extreme blue wing in our multi-configuration spectrum 
means that there is merely more CO emitting at higher 
velocities than at lower, which is indicative of a sheet like 
structure moving towards the observer. 

The line profiles of higher CO rot ational transitions for 
Bet elgeuse have been pub l ished inlKemper et al.1 ()2003l ) 
and lDe Beck et al.l (|2010t l. iDe Beck et al.1 (|2010T > present 
high spectral resolution (0.3125 MHz) line profiles for the 
CO(J = 2 - 1), (J = 3 - 2), and ( J = 4 - 3) transitions 
that were obtained with the James Clerk Maxwell Tele- 
scope (JCMT). For the CO(J = 2-1) transition the 
JCMT has an HPBW of - 20" and the profile appears 
similar to our multi-configuration profile over the same 
flux density extraction area (i.e., Figure [2]), with the ex- 
treme blue wing component being the dominant feature 
in both. This feature, which is emission from the S2 flow, 
becomes a less dominant component of the line profile at 
the higher CO( J =3-2) and CO(J = 4-3) transitions 
where the JCMT has an HPBW of - 13" and 8", respec- 
tively, and does not capture all of the S2 emission which 
is shown in Figure [3] to be extended at these velocities. 
Also, the higher rotational states (J « 10) will be popu- 
lated more by the higher excitation temperature (~ 200 
K) SI flow so these line profiles become dominated by 
emission from the slower SI flow. This is confirmed by 
our narrow SOFIA-GREAT CO(J = 12 - 11) line pro- 
file which is presented in the appendix of this paper, and 
also by a visual inspection of Herschel-HIFI archival line 
profiles of the CO(J = 6 - 5), (J = 10 - 9), and (J = 
16 — 15) transitions. 

The CO 4.6 /im ro-vibration li nes have been obser ved 
with the Phoenix spectrograph (jHinkle et al.1 fl998l ) by 
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iRvde et al.l (|1999t ) and lSmith et all (|2009D on the 2.1 m 

telescope at Kitt Peak and on the 8 m Gemini South 
telescope respectively. By assuming a Boltzmann pop- 
ulati on distribut i on fo r the ground rotational levels of 
CO. IRvde et~aTl (fl999h derived a mean excitation tem- 
perature of 38^5 K along the line of sight at a projected 
distance of 4" north of Betelgeuse. Our CARMA data 
suggests that the SI flow extends out to approximately 
this distance but Ryde et al.'s temperature is not in 
agreement with either of the linc-of-sight SI or S2 exci- 
tation temperatures of 200^0 K and 70±10K derived by 
iBernat et al.l (|1979f ). This discrepanc y may indica t e that 
the excitation is quite non- uniform. ISmith et all (|2009fi 
did not derive an excitation temperatures but used their 
4.6 (im spectra to reveal extended resonantly scattered 
CO emission out to ~ 2", a factor of two smaller than our 
SI radius. They observe emission over a velocity range 
of 30 km s _1 but two distinct flows are not detected. 
Mild (~ 20%) density inhomogeneities are reported but 
overall, their observations are consistent with an opti- 
cally thin and steady wind which is consistent with our 
findings. 



4.2. K 1 7699 A Spectra 

The S2 flow was first identifie d in high resolut i on KI 
and Nal absorption spectra by IGoldberg et al.l (|1975[ ) 
and subsequently re- observed multip le times over the 
next couple of years (|Goldberg||1979D . It is interesting 
to compare these and Bernat et al.'s (1979) CO line-of- 
sight absorption velocities with those from the CARMA 
emission spectra obtained at similar spectral resolutions 
and also to measure, the perhaps co-spatial line broad- 
ening of the K I S2 absorption feature. 

We have obtained KI 7698.98 A spectra using the 
cross-dispersed echelle spectrometers on the Harlan J. 
Smith 107 inch (2.7m) reflector at McDonald Obser- 
vatory. With 2 pixels per resolution element an R = 
A/AA = 200,000 and a R = 500,000 spectrum were 
obtained on 2007 March 25 and April 13, respectively. 
The spectra were wavelength calibrated with ThAr lamp 
lines and the lower resolution spectrum was checked by 
fitting six symmetric terres trial O2 lines in the same or - 
der using wavelengths from iBabcock fc Herzberd (|1948l ) . 
The 2 lines confirmed that the R = 200,000 cali- 
bration was good to better than 0.1 kms" . Upon 
cross-correlating the low- and high-resolution spectrum 
the high-resolution spectrum appeared redshifted by 
0.60 kms -1 , i.e., one resolution element, for which we 
do not have an explanation except to no t e that a similar 
offset has been reported bv lWeltv et al.l (|1994[ ). We use 
the cross-correlation to define the wavelength calibration 
of the R = 500, 000 spectrum and we adopt a systematic 
error of a sys = 0.2 kms" 1 . 

The high-resolution spectrum is shown in Figure [7] in 
the adopted stellar center-of-mass rest frame (14ad = 
+20.7 kms" ). The S2 feature is deep, well sepa- 
rated from the SI feature, and very well represented 
by an absorption model including the hyperfine split- 
ting from the abundant 39 K isotope. We adopt the 
KI 7698.9645 A line parameters compiled in iMortonl 



20 I ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 
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Figure 7. KI 7698.98 A profile (_R=500,000) for Betelgeuse ob- 
tained on 2007 April 13. The S2 flow outflow velocity is found to 
be 15.6 km s~ 1 (using a V ra d = +20.7 km s _1 ), slightly lower than 
IPlez & Lamberti 's lf2002T ) value of 18±2 km s" 1 . 

(|2003f Pl and find a heliocentric S2 absorption velocity of 
+5.15 ± 0.01(ler ran ) kms -1 with the uncertainty domi- 
nated the absolute wavelength scales. However, the most 
probable line-of-sight turbulent velocity is much more ac- 
curately determined with 0.62±0.02(ler ran )km s~ . From 
their 4.6 /im line profiles Bernat et al. (1979) found a 
very similar 1 km s" , these being the sum of thermal 
and non-thermal motions. These values place tight con- 
straints on the line-of-sight acceleration associated with 
the S2 flow. 

The K I spectrum also reveals a slight inflection in the 
observed line profile at +3.6 kms -1 (heliocentric) which 
may represent structure in the underlying photospheric 
profile or additional absorption in which case it has ~ 0.1 
the column density of S2 (A K i_s2 ^ 1-2 x 10 11 cm" 2 ). 
The S2 absorption velocity minima can be compared to 
those obtained by iGoldbergl (|1979l Figure 7) who mea- 
sured values between 1975 and 1978 of 4.2 ± 0.2 and 
5.0 ± 0.2 kms -1 and these differences may result from 
changes caused by radial velocity changes in the un- 
derlying photospheric spectrum. Bernat et al.'s 1979 
CO IR absorption observations reveal S2 heliocentric 
velocities of +4.94 + 0.30 kms" 1 (1979 March 6) and 
+4.60 ±0.04 kms" 1 (1979 April 14) with turbulent ve- 
locities of 4 and 1 kms" 1 for the SI and S2 features, 
respectively. 

In terms of the center-of-mass radial velocity of the 
star our KI feature implies an outflow velocity of 
+ 15.6 kms" 1 . The blue edge of our CARMA multi- 
configuration CO profile is estimated to be +15.4 km s" 1 
which suggests a dynamical associ ation with the CO S2 
flow and very close agreement with lBernat et all's (I1979T) 
CO ab sorption velocities listed above. IPlez &: Lambert! 
(|2002h have also estimated the radius and velocity of the 
suspected KI S2 flow using R = 110, 000 resolution long 

9 Note that this wavelength was recently revised and is 
0.44 kms -1 less that that adopted in the Goldberg studies. 
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slit spectra. They found a geometrically thin shell (1") 
with velocity of Vs2 = 18 + 2 kms" with a radius of 
55" which is much larger than the field of view of the 
CARMA spectra. Their long slit spectra show several 
smaller partial shells but it is not simple to directly asso- 
ciate the CO emission feature with one or more of these 
shells especially given the uncertainty in the ionization 
balances of CO and K I. It is possible that the 55" shell 
is associated with the inflexion caused by additional ab- 
sorption (and low column density) at a velocity slightly 
higher then S2 observed in our K I profile. 

5. CONCLUSIONS 

The two distin ct velocity components seen by 
IBernat eT al. (1979]) in CO absorption against the stellar 
spectrum at 4.6 /im have both been detected at 230 GHz 
for the first time. The first velocity component known as 

51 ha s an expansion velocity of 9 kms -1 (IBernat et al.l 
119791) and is detected in our high spectral resolution C 
configuration profile with the same blueshifted velocity 
(i.e., —9.0 kms -1 ) and with a larger redshifted outflow 
velocity of +10.6 kms -1 . The extended CARMA C con- 
figuration has a resolving out scale of ~ 6" and thus 
spatially filters almost all of the S2 emission leaving us 
with an approximate spectrum for the SI flow. An ex- 
treme blue wing of the CO spectrum appears in the D 
and E configuration spectra which we associate with the 

52 flow. The high spectral resolution multi-configuration 
spectrum is used to determine S2 outflow velocities of 
— 15.4 kms -1 and +13.2 kms -1 which is in good agree- 
ment with our K 1 7699 A line-of-sig ht S2 velocity and 
that reported bv IBernat et al.l (|1979( ). 

The low spectral resolution multi-configuration maps 
provide the first direct measurements on the spatial ex- 
tent of the S2 flow, which we derive to have a radius 
of 17"; a value that is higher than most previous esti- 
mates. We do not see a well defined outer edge for the 
SI flow but believe that it may extend out to a radius of 
~ 4". Previous single-dish observations of the CO line 
with small HPBWs do not show the classical resolved sig- 
nature of high emission at large absolute velocities and 
low emission at low absolute velocities for two main rea- 
sons. First, the SI flow is still unresolved in these single- 
dish observations and thus contributes emission and at 
the lower absolute velocities. As well as this, the multi- 
configuration CARMA maps show that the S2 emission 
is brighter in the higher absolute velocity maps than at 
lower absolute velocities and so when the emission from 
the fainter rings is neglected (i.e., when observed with 
a small HPBW) , the overall line profile does not change 
significantly. 

Assuming a mean outflow velocity of 14.3 km s _1 and 
9.8km s -1 for the S2 and SI flows respe ctively then their 
ages a re ~ 1100 yr and ~ 380 yr. Since iPlez k, Lambertl 
(|2002f) have detected KI out to 55" at a similar velocity 
to the CO S2 flow, then, assuming the CO and KI are 
coupled, there appears to be little or no further accel- 
eration in Betclgeuse's outflow once the S2 flow begins 
(which is somewhere greater than 4"). The composition 
and dynamics of the interface between SI and S2 remains 
unknown and future instruments such as the Atacama 
Large Millimeter/submillimctcr Array (ALMA) will pro- 
vide a greater understanding of this region. Higher spa- 



tial resolution, increased sensitivity, and excellent u — v 
coverage are needed to determine if the inner SI flow is 
discrete or just an extension of the current wind phase. 
Our SOFIA-GREAT spectrum shows that the higher ex- 
citation gas traces the slower SI component, and there- 
fore the high frequency bands of ALMA will preferen- 
tially trace the SI flow. Solutions to these remaining 
puzzles will broaden our knowledge of the evolutionary 
aspect of Betelgeuse's outflow and shed light into the 
driving mechanisms of M supergiant winds. 
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the states of California, Illinois, and Maryland, the 
James S. McDonnell Foundation, the Gordon and Betty 
Moore Foundation, the Kenneth T. and Eileen L. Norris 
Foundation, the University of Chicago, the Associates 
of the California Institute of Technology, and the Na- 
tional Science Foundation. Ongoing CARMA develop- 
ment and operations are supported by the National Sci- 
ence Foundation under a cooperative agreement, and by 
the CARMA partner universities. This work is based (in 
part) on observations made with the NASA/DLR Strato- 
spheric Observatory for Infrared Astronomy. SOFIA 
Science Mission Operations are conducted jointly by 
the Universities Space Research Association, Inc., under 
NASA contract NAS2-97001, and the Deutsches SOFIA 
Institut under DLR contract 50 OK 0901. GREAT is 
a development by the MPI fur Radioastronomie and 
the KOSMA/ Universitat zu Koln, in cooperation with 
the MPI fur Sonnensystemforschung and the DLR In- 
stitut fur Planetcnforschung. We thank Sarah Kennelly 
for providing us with archival Herschel spectra. This 
publication has emanated from research conducted with 
the financial support of Science Foundation Ireland un- 
der Grant Number SFI11/RFP.1 /AST/3064, and a grant 
from Trinity College Dublin. 

Facilities: CARMA, Smith and SOFIA. 

REFERENCES 



Airapetian, V. S., Ofman, L., Robinson, R. D., Carpenter, K., &; 

Davila, J. 2000, ApJ, 528, 965 
Altenhoff, W. J., Thum, O, & Wendker, H. J. 1994, A&A, 281, 

161 

Babcock, H. D., & Herzberg, L. 1948, ApJ, 108, 167 

Becklin, E. E., & Gehrz, R. D. 2009, in Astronomical Society of 
the Pacific Conference Series, Vol. 417, Submillimeter 
Astrophysics and Technology: a Symposium Honoring Thomas 
G. Phillips, ed. D. C. Lis, J. E. Vaillancourt, P. F. Goldsmith, 
T. A. Bell, N. Z. Scoville, & J. Zmuidzinas, 101 

Bernat, A. P., Hall, D. N. B., Hinkle, K. H., & Ridgway, S. T. 
1979, ApJ, 233, L135 

Bookbinder, J. A., Stencel, R. E., Drake, S. A., Simon, T., 
Linsky, J. L., & Florkowski, D. 1987, in Lecture Notes in 
Physics, Berlin Springer Verlag, Vol. 291, Cool Stars, Stellar 
Systems and the Sun, ed. J. L. Linsky & R. E. Stencel, 337 

Bracewell, R. N. 2000, The Fourier Transform and Its 
Applications, 3rd edn. (Colombus, OH: McGraw Hill) 

Carpenter, K. G., & Robinson, R. D. 1997, ApJ, 479, 970 

Danchi, W. C, Bester, M., Degiacomi, C. G., Greenhill, L. J., & 
Townes, C. H. 1994, AJ, 107, 1469 

Dc Beck, E., Decin, L., de Koter, A., Justtanont, K., Verhoelst, 
T., Kemper, F., & Menten, K. M. 2010, A&A, 523, A18 

Drake, S. A., Bookbinder, J. A., Florkowski, D. R., Linsky, J. L., 
Simon, T., & Stencel, R. E. 1992, in Astronomical Society of 
the Pacific Conference Series, Vol. 26, Cool Stars, Stellar 
Systems, and the Sun, ed. M. S. Giampapa & 
J. A. Bookbinder, 455 



CARMA CO(J = 2-1) Observations of Betelgeuse's CSE 



11 



Famaey, B., Jorisscn, A., Luri, X., Mayor, M., Udry, S., 

Dcjonghe, H., & Turon, C. 2005, A&A, 430, 165 
Glassgold, A. E., & Huggins, P. J. 1986, ApJ, 306, 605 
Goldberg, L. 1979, QJRAS, 20, 361 
— . 1984, PASP, 96, 366 

Goldberg, L., Ramsey, L., Testerman, L., & Carbon, D. 1975, 

ApJ, 199, 427 
Guesten, R. , et al. 2000, in Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Series, ed. 

R. K. Melugin & H.-P. Roser, Vol. 4014, 23-30 
Harper, G. M., Brown, A., & Guinan, E. F. 2008, AJ, 135, 1430 
Harper, G. M., Brown, A., & Lim, J. 2001, ApJ, 551, 1073 
Harper, G. M., Carpenter, K. G., Ryde, N., Smith, N., Brown, J., 

Brown, A., & Hinkle, K. H. 2009, in American Institute of 

Physics Conference Series, ed. E. Stempels, Vol. 1094, 868-871 
Hartmann, L., & Avrett, E. H. 1984, ApJ, 284, 238 
Hinkle, K. H., Cuberly, R. W., Gaughan, N. A., Hcynsscns, J. B., 

Joyce, R. R., Ridgway, S. T., Schmitt, P., & Simmons, J. E. 

1998, in Society of Photo-Optical Instrumentation Engineers 

(SPIE) Conference Series, ed. A. M. Fowler, Vol. 3354, 810-821 
Huggins, P. J. 1987, ApJ, 313, 400 

Huggins, P. J., Bachiller, R., Cox, P., & Forveille, T. 1994, ApJ, 
424, L127 

Kemper, F., Stark, R., Justtanont, K., de Koter, A., Ticlens, 
A. G. G. M., Waters, L. B. F. M., Cami, J., & Dijkstra, C. 
2003, A&A, 407, 609 

Knapp, G. R., & Morris, M. 1985, ApJ, 292, 640 



Knapp, G. R., Phillips, T. G., & Huggins, P. J. 1980, ApJ, 242, 
L25+ 

Lambert, D. L., & Vanden Bout, P. A. 1978, ApJ, 221, 854 
Lim, J., Carilli, C. L., White, S. M., Beasley, A. J., & Marson, 

R. G. 1998, Nature, 392, 575 
Morton, D. C. 2003, ApJS, 149, 205 
Plez, B., & Lambert, D. L. 2002, A&A, 386, 1009 
Rich, J. W., de Blok, W. J. G., Cornwell, T. J., Brinks, E., 

Walter, F., Bagetakos, I., & Kennicutt, Jr., R. C. 2008, AJ, 

136, 2897 

Ryde, N, Gustafsson, B., Hinkle, K. H., Eriksson, K., Lambert, 

D. L., & Olofsson, H. 1999, A&A, 347, L35 
Sanford, R. F. 1933, Contributions from the Mount Wilson 

Observatory / Carnegie Institution of Washington, 464, 1 
Scott, S. L., et al. 2004, in Astronomical Society of the Pacific 

Conference Series, Vol. 314, Astronomical Data Analysis 

Software and Systems (ADASS) XIII, ed. F. Ochsenbein, 

M. G. Allen, & D. Egret, 768-+ 
Smith, M. A., Patten, B. M., & Goldberg, L. 1989, AJ, 98, 2233 
Smith, N., Hinkle, K. H., & Ryde, N. 2009, AJ, 137, 3558 
Stebbins, J., & Huffer, C. M. 1931, Publications of the Washburn 

Observatory, 15, 178 
Taylor, G. B., Carilli, C. L., & Perley, R. A., eds. 1999, 

Astronomical Society of the Pacific Conference Series, Vol. 180, 

Synthesis Imaging in Radio Astronomy II 
Welty, D. E., Hobbs, L. M., & Kulkarni, V. P. 1994, ApJ, 436, 152 
Weymann, R. 1962, ApJ, 136, 844 



APPENDIX 



SOFIA-GREAT OBSERVATION 



As part of our larger multi-wavelength study of the CO surrounding Bctclgcusc we observed the star (PI: Harper ; 



ID 81_0005_1) with the German Receiver for Astronomy at Terahertz Frequencies ( GREA T: iGuesten et al.l 



2000) 



2009) 



instrument on NASA and DLR's Stratospheric Observatory for Infrared Astronomy (SOFIA; iBecklin fc Gehra 
2.5m airborne observatory. The 12 C 16 0( J = 12 — 11, 1.38 THz, 216.9 /im) line was observed to examine the dynamics 
of the higher excitation SI component. The observations were made during Flight 86 on 2011 November 10 at 13,100 m 
(~ 43,000 ft) when the star had an elevation of 45°. The HPBW was ~ 19" and the effective on source exposure time 
was 12 minutes. Due to technical difficulties during the flight the observations were obtained using a non-standard 
asymmetric chop sequence with a throw of 60". 

The scaled SOFIA-GREAT emission line profile binned to 1.3km s _1 is shown in Figure [8] along with the CARMA C 
configuration CO(J = 2 — 1) profile which spatially filters the S2 contribution. This figure shows that the J = 12 — 11 
profile reflects the slower moving SI flow with a width of ~ ±7.5 kms~ approximately centered on the stellar rest 
frame. This suggests that the J = 12 — 11 emitting plasma is not associated with the faster S2 component and is likely 
associated with the higher excitation SI plasma. Owing to uncertainties in the pointing accuracy during our GREAT 
observation we defer a discussion of the fluxes to a later time. 
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Figure 8. Scaled SOFIA-GREAT CO(J = 12 - 11) emission line profile binned to 1.3 kms -1 along with the unsealed CARMA C 
configuration CO(J = 2 — 1) profile. The star symbol at —2.5km s -1 marks the location of a weak terrestrial atmospheric emission feature 
which may contribute to the drop in flux at this velocity in the SOFIA-GREAT line profile. 



